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Abstract 

Application of an in-plane magnetic field to rippled graphene will make the system be a plane 
with randomly distributed vector potentials. Massless Dirac fermions carrying charges on graphene 
are scattered by the vector potentials and magnetoresistance is induced proportional to the square 
of amplitude of in-plane magnetic field Bd. Recently, Lundeberg and Folk proposed a formula 
showing dependence of the magnetoresistance on carrier density, in which the coefficient of B^ is 
given by a functional of the height-correlation function c(r) of ripples. In the present paper, we give 
exact and explicit expressions of the coefficient for the two cases such that c(r) is (i) exponential 
and (ii) Gaussian. The results are given using well-known special functions. Application of the 
present solutions in the vicinity of charge neutrality point should be careful because of the possible 
strong density inhomogeneity in rippled graphene. Here our analytic expressions are proposed, 
however, as trial interpolation formulas connecting the positive high-carrier-density regime and 
the negative one for magnetoresistance. Numerical fitting of our solutions to experimental data 
were performed. It is shown that the experimental data are well-described by the formula for the 
Gaussian height-correlation of ripples in the whole region of carrier density. The standard deviation 
Z of ripple height and the correlation length R of ripples are evaluated, which can be compared 
with direct experimental measurements. 
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I. INTRODUCTION AND MAIN RESULTS 



Graphene, a one- atom-thick allotrope of carbon [1, 2], has attracted much attention of 
theoretical physicists, since its electronic excitations are described by two-dimensional (2D) 
massless Dirac fermion (see, for instance, [3, 4]). Let v F be the Fermi velocity. Then, for a 
free fermion, the 2D massless Dirac equation is given by 



where % = y/—T, d x = d/dx, d y = d/dy, and ip{r) is a two-component wave function with 
the particle position r = (x,y). For a wave-number vector k = (k x ,k y ), the eigenvalues of 
energy are given by E = E± = ±hvpk, k = \k\, where the ± signs correspond to two values 
of the helicity of fermion h — ±1/2. The wave functions are determined as [3] 



Although the basic description of the system is very simple as shown above, physics 
of graphene is extremely rich, since the 2D Dirac equation is much modified by external 
electric and magnetic fields applied to the graphene plane as well as by changing geometric 
and topological features of the 2D surface [3, 4]. As a matter of fact, an ultrathin membrane 
is unstable in the 3D real space and in suspended graphene and in graphene placed on a 
substrate nanometer-scale corrugations, which are simply called ripples [5], are observed 
[6-12]. Since spatial distribution of such microscopic structures on the surface cannot be 
controlled, even when we apply uniform electric and/or magnetic fields, the rippled graphene 
will behave as a 2D system with randomly distributed scalar and vector potentials on it. 
Effect of ripples on electronic transport on graphene is an interesting subject both from 
theory and experiment [3, 4, 13-16]. 

Recently, magnetotransport measurements on rippled graphene have been reported, 
where a magnetic field is applied parallel to the averaged 2D plane of the graphene [17- 
20] (see also [21] for bilayer graphene). If the surface is a perfect plane, such an in-plane 
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where V is the volume of the system and 
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magnetic field B\\ does not affect any motion of 2D electron, since electronic motion couples 
only to the component of magnetic field normal to the 2D plane. But on real graphene, 
the magnetic field becomes to include a component normal to the surface depending on the 
local slope around ripples. Then the system has inhomogeneous out-of-plane magnetic fields. 
With respect to electromagnetic property, we will be able to regard the rippled graphene 
with B\\ as a system of traveling 2D massless Dirac fermions, which are randomly scattered 
by ripples [4, 17]. 

In the present paper, we study a formula of Lundeberg and Folk for the magnetoresistance 
Ap, which was derived from a model based on the Boltzmann transport theory [17, 18]. The 
relative height of ripple on the averaged 2D plane is denoted by h(r) for each 2D position 
r = (x,y). We characterize the height distribution by the two-point height correlation 

c(r;r ) = (h(r )h(r + r)), (4) 

which is defined as a statistical average of product of the height at r and that at the 
position separated from r by a displacement vector r. For simplicity, we assume that the 
distribution of h(r) is translation invariant and isotropic in this paper. By this assumptions, 
(4) becomes independent of r Q and it depends only on the distance of two points r = \r\; 
c(r; r ) = c(r). As a function of wave-number k = \k\ of fermion, we define 

poo 

g(k) = / rW(rk)c(r)dr, (5) 
Jo 

where 

W(z) = jf * Jo (2z sin sin 2 ^# (6) 

with the Bessel function Jq{z) (see Appendix A for special functions [22]). The signed 
carrier- density of graphene is denoted by n = n + — n_, which is the difference of the density 
of electrons n + and that of holes n_. We assume that the temperature is low and the state 
is degenerated so that the wave-number of fermion k is well-approximated by the Fermi 
wave- number kp. With counting contributions from the two Dirac points K and K' and 
spin, k ~ k F is related with n as (see, for instance, Sect. II. A in [3]), 



k = k{n) = a/ 7r\n\, (7) 

where the absolute value of n is used by the particle-hole symmetry; E± = ±hvpk. Let 
B\\ = \Bu\ and 6 be the angle between Bn and the direction of the current J along which 



the resistance p is measured. Set Ap = p(B\\) — p(0), the i?y -induced part of resistance, 
which we call magnetoresistance. Then the Lundeberg-Folk formula is given as 

nB 2 

Ap(n, 9, B\\) = ^f(sin 2 9 + 3 cos 2 9)g(k(n)). (8) 

Remark that this formula is not given in the published paper [17] of Lundeberg and Folk, 
but the corresponding equations are found in Appendix of its arXiv version [18]. Lundeberg 
and Folk introduced a correlation length R > for ripple-height correlation c(r) and argued 
the asymptotic behavior of g(k) in k <C 1/R and k ^> 1/R, and from the latter estimate the 
behavior of (8) in high carrier-density region was predicted as Ap ~ |n| -3 / 2 , \n\ 1. This 
|n|~ 3 / 2 -law was emphasized in their paper and has been used for analyzing experimental data 
[17, 19, 20], while it is an asymptotic law in |n| 3> 1 for a special case with the Gaussian 
height-correlation of ripples as clarified below in the present paper. Here we call the general 
equation (8) the Lundeberg-Folk formula. 

The main purpose of the present paper is to show that, for the two cases 

(i) c(r) = Z 2 e~ r / R (exponential), and 

(ii) c(r) = Z 2 e~ {r/R)2 (Gaussian), 

with a variance of height Z 2 = c(0) = (h(r ) 2 ) = (h 2 ), exact and explicit representations for 
the anisotropic magnetoresistance (8) are obtained. (For the Gaussian correlated disorder, 
see [23].) The results are the following. Let F(a, /3, 7; z) and F(a, 7; z) be the Gauss hyper- 
geometric function and the confluent hypergeometric function, respectively (see Appendix 
A for definitions and basic properties of special functions used in this paper). Then, for (i) 
the exponential height-correlation 

g(k) = AZKfF Q,|,2;-(2ti!) 2 

= ^[ A "' 2,:M »-TT(W £(2iM) ' (9) 

and for (ii) the Gaussian height-correlation 

9(k) = ^F( 3 -,2--(kR) 2 



7T{ZR) 



h{{kR) 2 /2) - h{{kRf/2) e-^ ' 2 , (10) 
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where K(z) (resp. E{z)) is the complete elliptic integral of the first (resp. of the second) 
kind, and I v {z) is the modified Bessel function of the first kind (see Appendix A). From 
these exact expressions, we can readily obtain the asymptotics of the magnetoresistance (8) 
in high carrier-density region as follows; For \n\ ^> 1, 



Ap(n,0,B\O 



sin 



( ZBw) 2 

2 9 + 3cos 2 9 y |n|~ 3/2 log(#|n| 1/2 ), for exponential (i) 



(11) 



I(sin 2 # + 3cos 2 £)^|^ 
4 hR 



n 



-3/2 



for Gaussian (ii). 



The asymptotic for the Gaussian case (ii) shown in the second line in (11) is exactly the 
same as Eq.(3) reported in [17]. On the other hand, there is a logarithmic correction to the 
|n|" 3 / 2 -law in the case (i) with the exponential height-correlation of ripples. 

For the low carrier-density region, we have found the following behavior for the magne- 



toresistance; when \n\ ~ 0, 



Ap(n,6,B\\ 



, , 9 A-kZRBw) 2 
(sin 2 9 + 3 cos 2 9) - — 

f v 

(ttZRBu) 2 



-(sin 2 fl + 3cos 2 #) 
4 



h 



9 

1 - -7ri? 2 |n| 
2 1 



-nR 2 \n 
4 



for exponential (i) 



for Gaussian (ii). 



(12) 



That is, p(n, 9, B\\) shows a cusp at the charge neutrality point |n| = and the dependence 
on n around this point is linear oc — \n\. In the vicinity of charge neutrality point, electronic 
transport will be affected by strong density inhomogeneity in rippled graphene, which is 
called electron-hole puddle [4]. Moreover, in low carrier-density region, \n\ 1, we have 
kR <^ 1 and the semi-classical Boltzmann transport theory is afraid to be broken-down. 
In order to give 'the 0-th order approximation', to which corrections by more elaborate 
approximations or full quantum theory will be added [4], however, here we try to apply our 
analytic results to the experimental data in the whole region of n. The result seems to be 
excellent if we assume the Gaussian height-correlation of ripples. See Figures 2 and 3 in 
Section IV. 

The paper is organized as follows. In Section II we give a derivation of the Lundeberg- 
Folk formula (8) from a semi-classical model using the Boltzmann transport theory. Section 
III is devoted to the proofs of our exact expressions (9) and (10) and their asymptotics (11) 
and (12). In Section IV, we discuss comparison of the present analytic results with the 
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experimental data measured by Wakabayashi and his coworkers. Appendix A is prepared 
for giving the mathematical formulas for special functions used in the text. In Appendix B 
an explanation is given for the linear relation between resistance and inverse of relaxation 
time based on the Boltzmann transport theory. 

II. DERIVATION OF LUNDEBERG-FOLK FORMULA 

As remarked just below Eq.(8), the general expression of Ap as a functional of the height- 
correlation function c(r) of ripples, (8) with (5) and (6), is found only in [18]. Although the 
essence of its derivation is given there, the description is very brief. For convenience, here 
we would like to complete the derivation of the Lundeberg-Folk formula (8) following [18] . 

We take the averaged 2D plane of graphene as the x-y plane. Let h(r) be the height of 
ripple at r — (x, y). 

A. Geometric vector potential 

Under the in-plane magnetic field in the x-direction, 



the geometric vector potential of Berry [24], A{r) = (A x (r), A y (r)), is calculated as [25] 



It gives a potential expressed by a 2 x 2 complex matrix for Dirac fermions with electric 
charge — e and the Fermi velocity vp, 

U(r) = —ev F A(r) ■ cr 
= ev F B\\h(r)a y 



B|| = (S||,0), 



(13) 



A x (r) = 0, 




B\\h(r). 



(14) 




(15) 



where cr = (a x , a y 



) are Pauli matrices. 
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B. Scattering matrix 



The scattering matrix for an interaction operator U is calculated by the Fermi golden 
rule as 

Sri PP Uk,k') = ^\{k'\U\k)\ 2 S(E(k') - E(k)) 

= 2 ^-\{k'\U\k)\ 2 8{k'-k), (16) 

where E(k) = hvpk. We note that 

(k'\U\k) = J d 2 r' J d 2 r (k'\r') (r'\U\r) (r\k) (17) 

with 

(r'\U\r) = U(r)5(r -r'), (18) 

where the potential U(r) is given by (15), S(r — r') = 5(x — x')5(y — y'), and (k'\r) = (r'\k'y. 
By using the two-component plane wave solution (2) for E = E + = hv^k, (17) becomes 

-ik' r .„,„.„,„{ I — i \ ] e^' 7 " I e~ i9k ^ 2 



(k'\U\k) = J d 2 r e -^={e ie *' l2 e- ie *' /2 )\ev F B\\h{r) 



iO V2V \ e ^/2 



= -^sm h (q) (19) 

with q = k' — k and 

h*(q) = J d 2 re- i( l r h(r). (20) 
The square of (20), \h(q)\ 2 = h(q)h*(q), is then calculated as 

\h(q)\ 2 = j d 2 r j d 2 r'e l( l< r '- r \(r)h(r') 

= V J d 2 re l( l- r ^ J d 2 r h(r )h(r + r), (21) 

where we have renamed the integration variable as r — > r and r' — r — > r. The quantity 
involving the second integral in the second line of (21), (1/V) / d 2 r h(r )h(r + r), is a 
volume average of a product h(r Q )h(r + r) with respect to r over the system. If the dis- 
tribution of the ripple height h(r) on graphene is translation invariant, this volume average 
will converge to the height correlation function (4) as V — > oo, by the self- averaging property 
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of usual random systems. Moreover, it can be written as c(r), if the system is also isotropic. 
Therefore, provided that V is large enough, (21) is equal to the Fourier transform of c(r) 



c(q) = J d 2 r e 1 ^ 7 'c(r), q = \q\ 



(22) 



multiplied by the volume V of the system; \h{q)\ 2 = Vc{q). Inserting the above results into 
(16), the volume factor V is cancelled out as desired, and we obtain the following expression 
of scattering matrix, 

2ne 2 v F B 2 



*Sripplc(&) k ) — 



sin 2 9k+ „ 9k ' c(q)S(k' -k), q = k' - k, q = \q\, (23) 



h 2 2 
which was given as Eq.(A.l) in [18]. 

By definition of 9 k for k = (k x ,k y ) given by (3) and its analogue Qy for kl = (k x ,k' y ), 
under the energy conservation k' = k, q x = k' x — k x = kicosOk' — cos 9k), q y = k' y — k y = 
k(sm0k> — sin^), and thus 



Q = \q\ = ^Qx + Qy 

= 2ky/l - COs(9 k - e w ) 



2k sin ■ 



Ok — Ok' 



(24) 



Then, when q = k' — k with the condition k' = k, (22) becomes 

roc i-2tt 

c(q) — I drr d(pc(r)t 
Jo Jo 



iqr cos <p 



OO /"27T 

dr rc(r) / dip exp 
o Jo 



2%kr sin cos ip 



(25) 



C. Boltzmann equation and scattering rates 

Lundeberg and Folk introduced a semi-classical model based on the Boltzmann transport 
theory for the system, in which the 2D electronic state with a wave-number vector k is 
described by probability distribution f(k). The time evolution of f(k) follows the Boltzmann 
equation (see, for example, [26, 27]) 



l- t + e -E.V k )f(k) 



df(k) 
dt 



(26) 



in the external electric field E, where we have assumed the spatial homogeneity of distribu- 
tion / and omitted the term v ■ V r f(k) in the LHS. The 'collision term' in the RHS, which 



describes the scattering of particles or holes by random vector potentials caused by ripples 
in the present system, is given by 

~df(ky 



dt 



= D[SJ](k) 



with the operator 



D[SJ](k) 



d 2 k' 



(27T) 



;S(k,k') f(k')-f(k) 



(27) 



(28) 



where S = S + hippie is the total scattering matrix. Here So is the zero field scattering 
matrix assumed to be in the isotropic form, S (k, k') = s (k, q)S(k — k'), and S^ppie is given 
by (23). We set S(k, k') = S(k, k')5(k - k') and rewrite (28) as 

k 



D[SJ](k) 



(2n)< 



d6 k ,S(k,k') f(6 k ,)-f(9 k ) 



(29) 



where f{9 k ) shows the dependence of f(k) on the angular component 9 k of k, when k = \k\ 
is given. In order to study ^-dependence of scattering induced by B\\, which is imposed in 
the x-direction, we introduce the following orthogonal bases, 



fxiPk) — — J=T, fy(8k) 



(30) 



0F ' JV ™ V^F ' 

satisfying 

/TV 
d0 k f a {0k)ffi{6k) = <W, a,(3 = x, y. (31) 

The -By-induced parts of scattering rates, which give the inverses of transport times enhanced 
by I? ||, are obtained as 

PIT 

(32) 



A VJW = - / d9 k f a (9 k )D[S ripple Je}(9 k ) 

for a, (3 = x, y, by the orthogonality relation (31) of the two bases (30). 

Now combining (32) and (29) with (23) and (25), we have multiple-integral expressions 
for Ar^ such that integrals are performed with respect to 9 k ,9 k /,r and ip. Change the 
integral variables as 9 k + 9 k > — > 9 + , 9 k — 9 k > — > 0_, and let 



L a ,p(k,r, 9. 



in 



d6 4 



2tt 



dip exp 



2ikr sin — cos a? 

2 Y 



9 + + 9. 



n 



9 + -9_ 



-n 



■ 2 ®+ 

sin — 

2 

9 + + 0_ 



a, (3 = x, y. Then 



Ar aA k "> = — 4^ff~ k J drrc ( r ) J d9_L at/3 (k,r,9_), a,p = x,y. 



(33) 



(34) 



By the Boltzmann transport theory (see Appendix B for an explanation), the resistances 
induced by By are related with as 

&PaA k ) = J^ Ar «S( fc )' <*,P = x, y- (35) 
We performed the integrals (33) and obtained the results 

L xx (k, r, 6L) = 3L yy (k, r, 6L) = -3n J (^2kr sin sin2 y, ( 36 ) 

I/ XJ ,(A;, r, 6L) = -^L^A;, r, 6L) = -n J (^2kr sin sin 6L, (37) 

where Jo(z) is the Bessel function with index v — (see Appendix A). 
Then (35) with (34) gives 

Ap xx (k) = 3Ap yy (k) = -^9(k), (38) 
Ap xy (k) = A Pyx (k) = 0, (39) 

where g(k) is given by (5) with (6). The 'threefold relation' between Ap xx and Ap yy is found 
in (38) [17, 18, 28]. 



D. Lundeberg-Folk formula 

Let E = (E x , E y ) be the in-plane electric field applied in order to measure resistance. If 
the current is written as J = (J x , J y ), then the resistance coefficients Ap a> p, a, j3 — x, y, are 
defined as 

E x \ _ I Ap xx Ap xy \ J x 
E v \ 

Ap yX Ap yy I \ Jy 

In the isotropic system, Ap xy = Ap yx = as in (39), then E x = Ap xx J x , E y = Ap yy J y . 
We have chosen the direction of B\\ in the x-direction as (13). Then, if we denote the 
angle between B\\ and J as 9, J x = J COS 6 , Jy — J sin# with J = \ J\ and the By-induced 
resistance Ap for the current J is given by 

E J 



(40) 



Ap = 



J 2 

^{Ap XX jl+A Pyy j 2 y) 



= Ap xx cos 2 9 + Ap, y y sin 2 9. (41) 

Then by the threefold relation (38) between Ap xx and Ap yy , the Lundeberg-Folk formula 
(8) is obtained. 
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III. PROOFS OF RESULTS 



A. Proof of Eq. (9) 



Let c(r) = Z 2 e r l R and use the formula (A3) for v — 0. Then (5) with (6) is written as 



„2\n 



n=0 



nl 



-M r (n)Afy(n) 



(42) 



with 



M r (n) 



drr 2n+l e -r/R = ^+2-^ + 3^ 







2?r 



2n+2 



c?0 sin 

2 



_ 2v^Fr(n + 3/2) 



(43) 



T(n + 2) 

We apply the duplication formula (A2) of Gamma functions and use the Pochhammer symbol 
(A5), then we have 



n=0 v ' 

2 f {(3/2) B } 2 (-(2tfi) 2 ) B 



\2\n 



ra=0 



(2)„ 



(44) 



which prove the first equality of (9) by the definition of the Gauss hypergeometric function 
(A6). 

From the recurrence relations (A8) and (A9), we obtain the following two equalities, 



f (|,3, 2; _ (2M)2 ) = _JL_ 



1; -(2*:B) 2 ) - F (i \,\;-{2kRf 



f(I,?, 1 ;-(2^) = tt ^ 

Combining them gives 



F[~l,l;-(2kR)A. 



(45) 



2(yfci?) 1 



F ( -, -, 1; —(2kR) 2 . 

1 2' 2' ' v W 1 + (2A;i?) 2 



1 



/■■(-i i l;-(2*i2) 2 



(46) 



By the formulas (A14), which express the complete elliptic integrals K(z) and E(z) by using 
the Gauss hypergeometric functions, the second equality of (9) is proved. 
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B. Proof of Eq.(lO) 



Let c(r) = Z 2 e ( r / R ) 2 and use the formula (A3) for v — 0. Then (5) with (6) is written 



as 



where 



,.2\n 



n=0 



n! 



(47) 



M r (n)= drr 2n+1 e- r2 ' R2 = -R 2n+2 n\ (48) 
jo 2 

and M^(n) is given as in (43). Following the similar procedure mentioned above, we have 

n(ZR) 2 v^(3/2) n (-(kRfr 



9(k) 



E 

n=0 



(2)„ n! 



(49) 



which proves the first equality of (10) by the definition of the confluent hypergeometric 
function (A7). 

From the recurrence relations (A10) and (All), we obtain the following equality, 

F(a + 1, r , z) = F(a, >y-l;z)+ V [ 1 F(a + 1, 7 + 1; z). (50) 
By using this equality for a — 1/2 and 7 = 2, we can rewrite the first line of (10) as 



<K*) = 



n(ZRf 



(51) 



By the formula (A15), which expresses the modified Bessel function I v {z) by using the 
confluent hypergeometric functions, the second equality of (10) is proved. 

C. Asymptotics in high carrier-density region 



Case (i): the exponential height- correlation 

In Gauss' integral representation (A12), if we change the variable as u = (w — 1)/ (4(kR) 2 ), 
we have 



Then we obtain the equality 
F( 3 -, 3 -,2;-(2kR) 2 



l+4(fcii) 2 / 1 

dw(w-l) 1 ' 2 (1 



4(kRY 



-1/2 



w^' 2 . (52) 



1 ( A 1 

4 + 



ir(kR) 3 \ (kR)< 



-1/2 



K 



(2kR)< 



1 + {2kRf 



E 



(2kRy 



1 + (2kR)< 



• (53) 
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By the asymptotics of the complete elliptic integrals [29] 

4 

K\z) ~ log 



z 



E(z) ~ 1, as z ->■ 1, (54) 

we have the behavior 

F (| | 2; -(2A;i?) 2 ) ~ -J^ [l g(4 A /l + (2kRf) - l] , as k -+ oo. (55) 

Through the relation (7) between k and n, the Lundeberg-Folk formula (8) with the present 
result (9) gives the first line of the RHS of (11). 

Case (ii): the Gaussian height- correlation 

In this case, we use the asymptotic expansion formula of the confluent hypergeometric 
function, 



n=0 



It gives 



r(«) ^ n! 

K^ 2; -H"^oW + Wi?n ' (57) 

and we obtain the second line of the RHS of (11). This result is exactly the same as Eq.(3) 
reported in [17]. 



D. Behavior in low carrier-density region 

Case (i): the exponential height- correlation 

By using the first two terms of the series (A6) defining the Gauss hypergeometric function, 
we obtain from (9) 

g(k) = n(ZR) 2 1 - j-(kR) 2 + 0((kR) 4 ) . (58) 
The first line of the RHS of (12) is concluded from (58). 
Case (ii): the Gaussian height- correlation 
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By using the first two terms of the series (A7) defining the confluent hypergeometric 
function, we obtain from (10) 



g(k) 



ir(ZRf 



1 - -(kR) 2 + 0{{kRf) 



(59) 



It gives the second line of the RHS of (12). 



IV. CONCLUDING REMARKS 



The Lundeberg-Folk formula (8) shows that the magnetoresistance induced by B\\ is 
proportional to B 2 , 

Ap(n,e,B\\)=a(n,e)Bl (60) 
The coefficient a(n, 9) is given by 

a(n, 9) = ^(sin 2 9 + 3 cos 2 6)g(y/Tr\n\), (61) 

in which the function g depends on the height-correlation function c(r) of ripples on 
graphene. In the present paper, we have clarified the dependence of Ap on c(r) for the 
two cases, (i) c(r) = Z 2 e~ r l R (exponential) and (ii) c(r) = Z 2 e~ <J / R ^ 2 (Gaussian), where 
Z 2 = c(0) = (h 2 ) is the variance of height of ripples and R indicates the correlation length 
of ripples. 

In the low carrier-density region, the so-called electron-hole puddles appear on the rippled 
graphene and the homogeneity assumption of system will be invalid. In order to handle 
the inhomogeneity of the carrier density around the Dirac point, the present Boltzmann 
approach should be corrected including higher- level approximations [4, 30, 31]. Away from 
the charge neutrality point, however, the semi-classical modeling based on the Boltzmann 
transport theory works well to describe the electronic transport properties of graphene. Here 
we consider our analytic expressions for Ap as 'trial interpolation formulas' connecting two 
high-carrier-density regimes, \n\ ^> 1; the positive- charge regime (n > 0) and the negative 
one (n < 0). 

Our expressions for Ap reported in this paper include two parameters Z and R. If 
experimental data are provided, we can evaluate the standard deviation of ripple height 
Z = a/ ih 2 ) and the correlation length of ripples R by numerical fitting of the data to our 
analytic expressions (9) or (10). 
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FIG. 1: The coefficient a(n,6) of the Lundeb erg- Folk formula (61) for 9 = with the present 
exact solution (9) for the exponential height-correlation of ripples is drawn. The parameters are 
chosen as Z = 0.531 [nm] and R = 8.51 [nm] by the best fitting to the experimental data given by 
Wakabayashi and his coworkers (indicated by dots). 



FIG. 2: The coefficient a(n,0) of the Lundeb erg- Folk formula (61) for 9 = with the present 
exact solution (10) for the Gaussian height-correlation of ripples is drawn. The parameters are 
chosen as Z = 0.376 [nm] and R = 17.0 [nm] by the best fitting to the experimental data given 
by Wakabayashi and his coworkers (indicated by dots). The curve fits the experimental data very 
well in the whole region of carrier-density n. 
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FIG. 3: Enlarged figure of Fig. 2 in the vicinity of the charge neutrality point (|n| = 0). Dots are 
the experimental data given by Wakabayashi and his coworkers. The bold curve showing a cusp at 
1 77, | = is the present analytic solution (10) for the Gaussian height-correlation with Z = 0.376 [nm] 
and R = 17.0 [nm]. In this low carrier-density region, the |n| _3/,2 -law shown by dotted curves is 
invalid, since it diverges as \n\ — > 0. In the very vicinity of the point \n\ = 0, the simple expression 
given in the second line of the RHS of (12) will work well, which is shown by thin black lines. 

Figures 1 and 2 show results of numerical fitting of the Lundeberg-Folk formula (8) with 
our exact solutions (9) and (10) for the exponential and the Gaussian height-correlations to 
the experimental data, respectively, when 9 = in the both cases. The experimental data 
are given by Wakabayashi and his coworkers [19, 20]. As we can see in these figures, the 
fitting in the high carrier-density region is more excellent in the Gaussian case (Fig. 2) than 
in the exponential case (Fig.l). 

For the vicinity of the charge neutrality point, |n| ~ 0, the fitting result for the Gaussian 
height-correlation is shown by Fig. 3, in which the |n| -3 / 2 -law for the high carrier-density 
region |n| ^> 1 (the second line of (11)) and the expression for |n| ~ (the second line of 
(12)) are also plotted for comparison. The obtained values of parameters by the present 
numerical fitting are Z = 0.376 [nm] and R = 17.0 [nm], respectively. These values should 
be compared with the direct experimental measurements [7, 12, 17, 19, 20]. In particular, it 
was reported in [20] the AFM measurements gave Z = 0.118 [nm] and R = 15.2 [nm]. See 
also experimental values cited in [17, 18] from [7, 12] and others. (We note that recent high- 
resolution measurements and analysis have revealed smaller-sized roughness of the substrate 
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and graphene surfaces [32].) 

We hope that the our expressions will be used for systematic analysis of experimental 
data also for different values of 9. In the present paper, we have assumed translation 
invariance and isotropy for distribution of ripple height. Even in the high carrier-density 
regions, apart from the effect of electron-hole puddles, further calculation for general setting 
of ripple-height-correlations will be an interesting future problem. 
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Appendix A: Definitions and formulas of special functions 

Here we list the definitions and formulas of special functions used in the text (see, for 
instance, [22]). The gamma function T(z) is defined by 



When n e N = {0, 1, 2, . . . }, T(n + 1) = n\ (we set 0! = 1). The following equality is called 
the duplication formula, 



for a complex number z which is not negative real. Then the modified Bessel function of 
the first kind is defined as 




(Al) 




(A2) 




(A3) 




:/2 J u {ze wi/2 ) (-7i < aigz < tt/2) 
l2 J v {ze-^ 2 ) (n/2 < argz < vr) 



(A4) 
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For n G No, the Pochhammer symbol is defined as 

(a) n = a(a + l)(a + 2) • • • (a + n - 1), n G N = {1, 2, 3, . . . }, 

(a)o = 1. (A5) 



The Gauss hypergeometric function is defined by 

(a 

{l)n n\ ' 



F {a ,p, r ,x) = £®gk* (A6) 



n=0 

and the confluent hypergeometric function is defined as 



oo 



F(a,r,z) = Jam F{a,p,y,z/p) 

' ' o 

(a) n z 

n=0 

The following recurrence relations for F(a,j3, , y; z) are known, 



Vf^T (A7) 
^ (7)n n\ 



7[F(a, + 1, 7; *) - F(a, p, 7; 2)] = c^F(a + 1, p + 1, 7 + 1; 2), (A8) 
a(l - z)F(a + 1, /3, 7; *) + [7 - 2a + (a - p)z]F(a, p, 7; z) 

+(a->y)F(a-l,P,r,z) = 0. (A9) 

For F(a,p;z), 

zF(a + 1, 7 + 1; 2) = 7 [F(a + 1, 7; z) - F(a, 7; 2)], (A10) 

aF{a + 1, 7 + 1; z) = (a - i)F(a, 7 + 1; z) + ~fF(a, 7; z) (All) 

are satisfied. The Gauss hypergeometric function has the following integral expression, 

F(a, P, 7; z) = T{ ^_ /3) f Q du vP-\l - uy-f>-\\ - uz)-«. (A12) 

The complete elliptic integrals of the first kind, K(z), and of the second kind, E(z), are 
defined as 

K(z) = f 72 de = f 1 du 

1 ' ~ Jo Vl-z^m 2 e Jo v/(l-« 2 )(l-^ 2 )' 
E{z) = / VI- z 2 sin 2 6d6, (A13) 



respectively. They are expressed by using the Gauss hypergeometric functions as 

^) = fF(i,I,l;/), BW = ^ (-!,!,!;/). (AM) 
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The modified Bessel function (A4) is expressed by using the confluent hypergeometric func- 
tions as 

I v {z) = Y^yfF (u + \ 2v + 1; -2s) . (A15) 
Appendix B: Linear relation between p and r 1 

Let fo(k) = [e^ E ^~^^ kBT ^ + be the equilibrium Fermi distribution function at tem- 
perature T and chemical potential /i. The relaxation time r{k) is phenomenologically intro- 
duced as 



df(k) 
dt 



with 5f(k) /(/.-)• ./„(/.•). (Bl) 



r{k) 



When we consider the long-term limit, df(k)/dt-term in the LHS is zero, and we have 

Sf(k) = - e - T {k)E-V k f{k) 

~ -2_ T ( k )E-V k fo(k). (B2) 

In low temperature limit, E = hv F k and E ■ Vkfo(k) ~ hv F \E\df (k)/dE ~ |.E|5(A; — &f), 
where k F denotes the Fermi wave-number. Then in the 2D plane the current is calculated 
as follows, provided J and E are parallel; 

J=\J\= J ^- 2 (-ev F )5f(k) 

—k(-ev F ) [--r(k)\E\5(k - k F )) 

= ±^IlM\ E \. (B3) 
2'nh 

Since resistance p is defined as J = p^\E\, we have the linear relation between p and r _1 
as p = ct^ 1 with the coefficient c = 2nh/ \e 2 v F k F ). In Eq.(35) k F is written as k. 
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